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^Nj [ Context. The study of the signatures of multimodality in groups and clusters of galaxies, an environment for most of the galaxies in 

the Universe, gives us information about the dynamical state of clusters and about merging processes, which affect the formation and 
evolution of galaxies, groups and clusters, and larger structures - superclusters of galaxies and the whole cosmic web. 
Aims. We search for the presence of substructure, a non-Gaussian, asymmetrical velocity distribution of galaxies, and large peculiar 
r \ ' velocities of the main galaxies in clusters with at least 50 member galaxies, drawn from the SDSS DR8. 

, Methods. We employ a number of 3D, 2D, and ID tests to analyse the distribution of galaxies in clusters: 3D normal mixture mod- 

^ ■ elling, the Dressler-Shectman test, the Anderson-Darling and Shapiro- Wilk tests, as well as the Anscombe-Glynn and the D' Agostino 

CLf tests. We find the peculiar velocities of the main galaxies, and use principal component analysis to characterise our results. 

Results. More than 80% of the clusters in our sample have substructure according to 3D normal mixture modelling, and the Dressler- 
Shectman (DS) test shows substructure in about 70% of the clusters. The median value of the peculiar velocities of the main galaxies 
in clusters is 206 km/s (41% of the rms velocity). The velocities of galaxies in more than 20% of the clusters show significant non- 
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ABSTRACT 



o 



^ ■ Gaussianity. While multidimensional normal mixture modelling is more sensitive than the DS test in resolving substructure in the sky 

^ ' distribution of cluster galaxies, the DS test determines better substructure expressed as tails in the velocity distribution of galaxies 

(possible line-of-sight mergers). Richer, larger, and more luminous clusters have larger amount of substructure and larger (compared 
to the rms velocity) peculiar velocities of the main galaxies. Principal component analysis of both the substructure indicators and the 
^ ■ physical parameters of clusters shows that galaxy clusters are complicated objects, the properties of which cannot be explained with 

' a small number of parameters or delimited by one single test. 

, Conclusions. The presence of substructure, the non-Gaussian velocity distributions, as well as the large peculiar velocities of the 

^>«^ . main galaxies, shows that most of the clusters in our sample are dynamically young. 

Key words. Cosmology: large-scale structure of the Universe; galaxies: clusters: general 
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^ ' 1. Introduction and evolution of galaxies dBird & BeersI [1991 iPinknev et al] 



1996; Knebe & Miiller 2000). These mergers have shaped the 



^ Most galaxies m the Universe are located m groups and properties of galaxies in groups and clusters, as, e.g., the well 

• rH . clusters of galaxies, which themselves reside m larger sys- known morpho l ogical segregation effect (Einastoetll 

p><; ; tems - m superclusters of galaxies or m filaments crossing [Dresslei' W; 'Einasto & Einastd ^1987'; Berrier^ 
H ■ under d ense regions betwee n sup er clusters (iJoeveer et alJ 



1974 



2009; 



— --n — 7 ". — I — I — [Huertas-Companv et al. 2009). Substructure affects estimates of 

1978^, Gregorv&ThompsonI [19781; [Zeldovich et al.[ [1982|; several cluster characteristics , the dynamic al mas s and mass 



^apparent etal. 1986). In the ACDM concordance cos- j^.j^gj^^ ^^^^^ ^^^^^ ^j^ers dBiviano et al.l 12005 iNiemi et al.) 

mological model the sti-uctures forming the cosmic web 2007; Piff^aretti & Valdai-nini 2008; Holopainen etalj [2008[; 

grow bi^ierarc^cal clustering driven by gravity (see, e.g., , white et ah ,2010c ,Power et ah ,201 L and references therein). 

[Loebl [20021 [20081 and references therein). The present-day detailed knowledge of the properties of clusters of galaxies is 

dynamical state of clusters of galaxies depends on their ^^^^^^ comparison of observations with N-body models 
formation history. Signatures of multimodality in the distri- formation and e volution of cos mic str uctures, to test the 

bution of galaxies in clusters (the presence of substructure, cosmological models (iThomas et alJ [T998L lArava-Melo et all 

several galaxy associations within clusters, non-Gaussian 200% 
velocity distributions of galaxies, and large peculiar velocities 

of the main galaxies) are indicators of former or ongoing To seai-ch for signatures of multimodality in galaxy clus- 

mergers in groups and clusters, which affect the formation ters a number of 3D, 2D, and ID method s have been proposed 

including the Dressler-She ctman test dDressler & ShectmanI 

Send offprint requests to: M. Einasto 119881: iKnebe & Miille3l2OO0l) . the hierachical clustering method 
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Sema& Gerba| Il996t iDurret et al] l2010h . wavelet analysis 
Flin & Krvwult"2006'), multidimensional normal mixture mod- 
elling ( Fraley & Raftery 2006), see also Einasto et al. (2010), 
and many others (' see iPinknev et al.l [1996, for a review). 
iPinknev et al.l d 19961) showed that it is preferable to use several 
methods to search for substructure in clusters, since their sensi- 
tivity to different signatures of multimodality is different. 

Several studies have shown the presence of sub- 
structure in poor and rich gr oups and clusters of 



_galaxies (ISola nes et all '1999"; 'Oeaerle & Hilll 1200 U 
Kolokoti-onis et al.l 120011; FBursett et al. 2004; Boschi n et al 
2006; Fhn & Ki-ywult 2006; Bai-rena et al. 2007; Hwang & Led 



. 2007 - Boschin et al. 2008.; Hou et al.. 2 009; Ven nik & Hopd 
2OO9I; lAguerri & Sanchez-JanssenI l2010ri Pimbblet et all 1201 U 



Ribeiro et al.l 201 1 : iMartmez & Zandivarez 2012; 'Hou et al 



2012), and in X-ray clu sters (Ramella et al. 2007; Owers et al. 
2009blai ; iBohringer et all I2010() . ITovmassian & PlionisI (120091) 



studied the properties of poor groups from the SDSS survey 
and showed that many groups of galaxies are not presently in 
a dynamical equilibrium, but at various stages of virialization. 
[Niemi et al. (2007) showed that a significant fraction of nearby 
groups of galaxies are not even gravitationally bound systems. 

In this paper we study the multimodality in rich clusters 
drawn from the SDSS DR8. We use data about 109 clusters with 
at least 50 member galaxies, this is one of the largest samples 
of rich clusters analysed for substructure so far. We employ a 
number of 3D, 2D, and ID methods to analyse the distribution 
of galaxies in clusters. With the principal component analysis 
we characterise the results of different tests simultaneously, and 
study the relations between the multimodality of clusters and 
their physical properties. We present lists of unimodal and multi- 
modal clusters. In Sect.|2]we describe the data we used. In Sect. [3] 
we describe the methods to search for signatures of multimodal- 
ity, and we apply them in Sect. |4]to study the properties of clus- 
ters. We discuss the results and draw conclusions in Sect.|5] 

We assume the standard cosmological parameters: the 
Hubble parameter Hq - 100 h km s ' Mpc the matter den- 
sity Qn, - 0.27, and the dark energy density Qa - 0.73. 



2. Data 

We used the MAIN galaxy sample of the 8th data release of 
the Sloan Digital Sky Survey (Aiharaet al. 2011) with the ap- 
parent r magnitudes r < 17.77, and the redshifts 0.009 < 
z < 0.200, in total 576493 galaxies. We corrected the redshifts 
of galaxies for the motio n relative to the CMB a nd computed 
the co-moving distances (' Martinez & Saarl |2002|) of galaxies. 
The absolute magnitudes of galaxies were determined in the 
r-band (Mr) with the fc-corrections for the SDSS galaxies, cal- 
culated using the KCOR RECT algorithm dBlanton et al.ll2003at 
iBlanton & Roweisll2007h . In addition, we applied evolution cor- 
rections , using the luminosity evolution model of Blant on et al.l 
(I2003bh . The magnitudes correspond to the rest-frame at the red- 
shift z - 0- More deta ils on similar da ta reduction for the SDSS 
DR7 can be found in lTago et al.l (12010 . hereafter TIO). 

We determine groups of galaxies using the Friends-of- 
Friends (Fo F) cluster ana lysis method introduced in cos- 



mology by iTurner&GottI (Il976h : IZeldovich et al (Il982h : 

inP'a 



lHuchra& Gelled (1 1982h . and modified in^ago et al.1 liom and 



in TIO. A galaxy belongs to a group of galaxies if this galaxy has 
at least one group member galaxy closer than a linking length. 
In a flux-limited sample the density of galaxies slowly decreases 
with distance. To take this selection effect into account properly 



when constructing a group catalogue from a flux-limited sam- 
ple, we rescaled the linking length with distance, calibrating the 
scaling relation by observed groups (see TIO for details). As a 
result, the maximum sizes in the sky projection and the velocity 
dispersions of our groups are similar at all distances. This shows 
that distance-dependent selection effects have been properly ac- 
counted for. Our catalogue contains 77858 groups with at least 
2 member galaxies. The richest groups in our catalogue corre- 
spond to rich clusters of galaxies. The details and availablility 
of the group catalogue based on SDSS DR8 are described in 
ITempel etatl (2012). 

In flux-limited samples galaxies outside the observational 
window remain unobserved. To calculate the total luminosities 
of groups we have to take into account the luminosities of these 
galaxies as well. For that, we multiply the observed galaxy lu- 
minosities by the luminosity weight Wd- The distance-dependent 
weight factor Wj was calculated as follows: 



L«(L)dL 
j';Ln(L)dL 



(1) 



where Li,2 - LqIO'^ '*^*^©"*''*' are the luminosity limits of the 
observational window at a distance d, corresponding to the ab- 
solute magnitude limits of the survey Mi and M2; we took 
Mq - 4.64 mag in the r-band (Blanton & Roweis 2007), and 
n(L) is the galaxy luminosity function. Owing to their peculiar 
velocities, the distances of galaxies are somewhat uncertain; if 
the galaxy belongs to a group, we used the group distance (mean 
distance of galaxies in a group) to determine the weight factor. 
Detailed description, how the weig ht factor and luminos ity func- 
tion are calculated can be found in ITempel et al.l (1201 ll) . 

In the group catalogue the main galaxy of a group is defined 
as the most luminous galaxy in the r-band. We use this definition 
also in the present paper 

Next we select clusters for this study. The larger 
the number of galaxies in clusters, the more reliable is 
the analysis of substructure a nd th eir velocity distribu- 
tion. Aguerri & Sanchez-Janssen (1201 Oh select clusters with at 
least 30 member galaxies for substructure study. However, 
iBoschin et afl (12008 ) argue that clusters with about 30 mem- 
ber galaxies are too small for the analysis of substructure (see 
also the discussion about the substructure statistics in small sam- 
ples in Biviano et al. 2006), although clusters with smaller num- 
bers of galaxies have been studied for substructure (see, for 
example. lSoianes et al.l[T999t iBoschin et afll2008t iRibeiro et all 
1201 ll) . Another problem arises with the selection effects: in the 
group catalogue the richness of groups decreases rapidly at dis- 
tances D > 340 Mpc owing to the use of a flux-limited sam- 
ple of galaxies (TIO). At distances smaller than 120 /i"' Mpc the 
sample includes nearby exceptionally rich clusters which cor- 
respond to well-known Abell clusters (the Coma cluster, rich 
clusters i n the Herc ules supercluster and others, see TIO and 
lEinasto et al.ll201 Ibl) . These clusters have to be analysed sepa- 
rately. Therefore we chose for the present analysis clusters with 
at least 50 member galaxies in the distance interval 120 Mpc 
< D < 340 h-^ Mpc (redshift interval 0.04 < z < 0.12). This 
sample includes all clusters from the SDSS DR8 with at least 50 
member galaxies from our catalogue in this distance interval, in 
total 109 clusters. Figure [T] shows the richness of clusters in our 
sample vs. their distance. 

We cross-identify groups with Abell clusters, which have ac- 
quired a role of a reference system for rich clusters. Contrary to 
expectations, cross- identification of (rich) SDSS DR8 groups 
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with Abell clusters is not straightforward. Problems arise both 
due to different group/cluster finding procedures as well as dif- 
ferent algorithms for components/subgroups. The Abell clus- 
ters have a constant linear radius (1.5 Mpc), while the DR8 
groups obtained by a FoF procedure have various linear sizes. 
An Abell cluster may consist of several subclusters and/or may 
be the result of projections of groups from different distance (e.g. 
[Pimbblet et al. 201 1). These facts make cross-identification dif- 
ficult. We consider a group identified with an Abell cluster, if 
the distance between their centres is smaller than at least the lin- 
ear radius of one of the clusters, and the distance between their 
centres in the radial (line-of-sight) direction is less than 600 km/s 
(an empirical value). As a result one group can be identified with 
more than one Abell clusters and vice versa. This can be seen in 
online tables, where we give data on the clusters, and the results 
of the tests (Table [1] and Table |2|. In Fig. [2] we show the dis- 
tributions of cluster redshifts and total luminosities, as given in 
Table □ 
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Fig. 1. Richness of clusters vs. their distance. 
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Fig. 2. Distribution of cluster redshifts (upper panel) and total 
luminosities (lower panel). 



3. Methods 

In this section we decribe the methods applied in this paper to 
analyse the multimodality of galaxy clusters. 



3. 1 . Multidimensional normal mixture modelling with Mclust 

To search for possible components in clusters, we em- 
ploy multidimensional normal mixture modelling based on 
the analysis of a finite mixture of distributions, in which 
each mixture component is taken to correspond to a dif- 
ferent group, cluster or subpopulation. The most common 
component distribution considered in model-based cluster- 
ing is a multivariate Gaussian (or normal) distribution. To 
model the collection of components, we apply the Mclust 
packa ge for classification and clustering dFralev & Raftervl 
2006") from R, an open-source free statistical environment 
developed under the GNU GPL (Ihaka & Gentleman 199^ 
http://www.r-project.org). This package searches for an 
optimal model for the clustering of the data among models with 
varying shape, orientation and volume, finds the optimal number 
of components, and the corresponding classification (the mem- 
bership of each component). Mclust calculates for every galaxy 
the probabilities to belong to any of the components. The un- 
certainty of classification is defined as one minus the highest 
probability of a galaxy to belong to a component. The mean un- 
certainty for the full sample is used as a statistical estimate of 
the reliability of the results. 

We tested how the possible errors in the line-of-sight posi- 
tions of galaxies affect the results of Mclust, shifting randomly 
the peculiar velocities of galaxies 1000 times and searching each 
time for the components with Mclust. The random shifts were 
chosen from a Gaussian distribution with the dispersion equal to 
the sample velocity dispersion of galaxies in a cluster The num- 
ber of the components found by Mclust remained unchanged, 
demonstrating that the results of Mclust are not sensitive to such 
errors. 

We also performed substructure analysis with Mclust using 
only the sky coordinates of the clusters as an additional 2D test. 



3.2. DS test 

Another diagnostic for su bstructure in clusters is the Dressler- 
Shectman (DS or A) test dPressler & Shectman|[T988l) . The DS 
test searches for deviations of the local velocity mean and dis- 
persion from the cluster mean values. The algorithm starts by 
calculating the mean velocity (v'locai) and the velocity dispersion 
(ciocai) for each galaxy of the cluster, using its n nearest neigh- 
bours. These values of local kinematics are compared with the 
mean velocity (v^) and the velocity dispersion (cTc) determined 
for the entire cluster of Ngai galaxies. The differences between 
the local and global kinematics are quantified by 

= (« -I- l)/crl [(Viocal - Vcf + (o-iocal - (Tcf] . 

The cumulative deviation A - 25, is used as a statistic for 
quantifying (the significance of) the substructure. The results of 
the DS-test depend on the number of local galaxies n. Here, we 
ha ve studied the substructure by using n - y/Niyai, as suggested 
bv lPinknev etal] ([1996). 

If the cluster velocity distribution is close to a Gaussian, then 
A will be of order Ngai. For the case of non-Gaussian velocities, 
A can differ significantly from Nm i, even if there is no subclus- 
tering dDressler & Shectmanll988h . Therefore, the A statistic for 
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each cluster should be calibrated by Monte Carlo simulations. 
In Monte Carlo models the velocities of galaxies are randomly 
shuffled among the positions, which effectively destroys any true 
correlation between the velocities and positions. We ran 25000 
models for each cluster and calculated every time Asim. The sig- 
nificance of having substructure (the /j-value) can be quantified 
by the ratio A^(Asim > \,bs)/Nsim - the ratio of the number of 
simulations in which the value of A is larger than the observed 
value, and the total number of simulations. The smaller the p- 
value, the larger is the probability of substructure. 

3.3. a test 



With the a test developed by IWest & BothunI (119901) we anal- 
yse correlations between positions and velocities of galaxies and 
search for a region which differs from the overall distribution. 
This test gives a measure of the centroid shift for the cluster 
galaxies. The idea of the test is to search for a region that shows 
correlation between the positions and velocities that differs from 
that for the overall galaxy distribution. The method consists of 
a five different steps. At first the centroid of the system (xcjc) 
is calculated as a global mean value. Next each galaxy / gets 
a weight w, - l/V,, where V, is the line-of-sight rms velocity 
calculated using the velocities of the galaxy and its 10 nearest 
members. For each galaxy 10 nearest neighbours are chosen in 
the velocity space and new centroid (x^, y'j is calculated using 
the weights: 



■t-.ll v^ll 

■t-.li ' -^c ~ „ii 
2,=i Wi Y.i=x Wi 



(2) 



In the fourth step it is calculated how much this new centroid 
differs from the global centroid: 



a,- 



(xc - + (jc - 



Finally the average value of a, for all galaxies is defined as the a- 
value. This is a measure of how much the centroid of all galaxies 
shifts as a result of local correlations between the positions and 
velocities of galaxies. 

3.4. p test 

We also employ 2D tests which use information about the sky 
positions of galaxies in clusters: the test which stu dies the 
asymmetry in the galaxy distribution (IWest et al.lll988h and 2D 

normal mixture modelling with Mclust. 

The p test, presented in IWest et all (Il988h . is a test for the 
asymmetry in the galaxy distribution. At first, the mean distance 
di for each galaxy / to its five nearest neighbours is calculated. 
Then a point diametrically opposite to this point is chosen and 
the mean distance do in this neighbourhood is calculated as be- 
fore. The asymmetry is then measured by the /3-value: 



P = logiQ(do/di) 



(4) 



To evaluate the asymmetry the average value < yS > over 
all galaxies is calculated and any deviation from < /? 
signifies asymmetry and possible substructure. According to 
iPinknev et al.l (1 19961) the p test is sensitive to the mirror asym- 
metry, but not to the deviations from the radial symmetry. The 
p-values for the a and p tests are calculated with Monte Carlo 
simulations as for the DS test using 25000 random sampUng of 
galaxy coordinates in clusters. 



3.5. 1D tests 

One indicator of multimodality in clusters is the deviation of the 
distribution of the galaxy velocities in clusters from a Gaussian. 
To analyse the distribution of the velocities of galaxies we use 
several ID tests. We tested the hypothesis about the Gaussian 
distribution of the peculiar velocities of galax i es in clusters 
with the Shapiro- Wilk normality test dShapirol Il965h . which 
is considered the best for small samples. We also employed 
the Ander s on-D arling test, which is very reliable according to 
iHou et af] (l2009b . We calculated the kurtosis and the skewness 
of the peculiar velocity distributions, and used these to test for 
the asymmetry of the distributions of galax y velocities. We used 
the Anscombe-Glynn test for the kurtosis (Anscombe & Glvni 
1983), and the D'Agostino test for the skewness (D'Agostiri 
197()j), from the R package moments by L. Komsta and F. 
Novomestky. 

There is, of course, no reason to assume that the velocity dis- 
tributions in galaxy clusters should be exactly Gaussian. In fact, 
kinematic models show that the shape of the velocity distribution 
is defi ned by the ratio of different types of galaxy orbits (iMerrittI 
'1987"). The tests described above serve mostly to check if the 
velocity distribution is unimodal and symmetrical. Historically, 
they have been selected because of the easy availability of statis- 
tical tests for the Gaussian distribution; we use them to be able 
to compare our results with those obtained earlier 

In virialized clusters galaxies follow the cluster potential 
well. If so, we would expect that the main galaxies in clusters lie 
at the cent res of groups (group haloes) and have sniall peculiar 
veloc ities (lOstriker & Tremaineiri975l: iMerrittI [19841 iMalumuthI 
Il992h . Therefore the peculiar velocity of the main galaxies in 
clusters is also an indication of the dynamical state of the cluster 
(Coziol et al. 2009). We calculate the peculiar velocities of the 
main galaxies, Ivped, the normalised peculiar velocities of the 



main galaxies, Vp 



(3) main galaxies in the clusters and subclusters. 



/cr,,, and analyse the location of the 



3.6. Principal component analysis 

We employ the principal component analysis (PCA) to analyse 
the results of all tests simultaneously, and to study the relations 
between the multimodality indicators and the physical param- 
eters of clusters. The aim of the PCA is to study the relations 
between the parameters and, if possible, to find a small num- 
ber of linear combinations of correlated parameters to describe 
most of the variation in the dataset with a small number of new 
uncorrelated parameters. The PCA transforms the data to a new 
coordinate system, where the greatest variance by any projec- 
tion of the data lies along the first coordinate (the first principal 
component), the second greatest variance - along the second co- 
ordinate, and so on. There are as many principal components as 
there are parameters, but often only the first few are needed to 
explain most of the total variation. 

The principal components PC/ (/ e N, / < Mot) are linear 
combinations of the original parameters: 



' 'lot 

PC/ = ^fl(^);yi, 



(5) 



k=l 



where -1 < fl(fc), < 1 are the coefficients of the linear trans- 
formation, Vk are the original parameters and A^tot is the num- 
ber of the original parameters. In the analysis the parameters 
are standardised - they are centred on their means, Vk - Vk, 
and normalised, divided by their standard deviations, cr(y^). The 



4 



M. Einasto et al.: Substructure of clusters 



PCA is suitable tool for exploratory data analysis, to study si- 
multaneously correlations between a large number of parameters 
(see Einastoetjil. 201 la^ for references about applications of the 
PCA in astronomy). 



4. Results 

We show in Table [T] the numbers and fractions of clusters with 
substructure, and with the p-values showing statistically signif- 
icant deviations from Gaussianity in their galaxy distributions 
ip < 0.05). 



Co mparison with other studies shows that iRamella et al.l 
(I2007h found substructure in 73% of X-ray clusters which 
is very cl ose to what we fo u nd by 2D nor mal mixture 
modelling. ISolanes et alJ (Il999h lOegerle & Hilll and 



lAguerri & Sanchez-JanssenI (1201 Ol) found with the DS test that 
30-50 % of clusters have significant substructure. lFlin & Krvwul3 
(1200 6') detected substructure in about 1/3 of Abell clusters stud- 
ied by them using wavelet analysis. About 20-30% of clusters 
have non-Gaussian ve locity distributions, as shown with the AD 
test or other ID t ests ("Solanes et al."l999';'Oe2erle & Hill"200jt 
Hou et alJ 120091; [Rib eiro et al. 2011; Martinez & ZandivareS 
20T2I)! 



Table 1. Results of the tests. 



Population 




% 




% 




% 




All 




A'compSD 


= 1 


^ comp3D 


> 1 


All clusters 


109 




17 


16 


92 


84 


^comp2D > 1 


82 


75 


3 


18 


79 


86 


Vpec > 250 


48 


44 


4 


24 


44 


48 


Vpecr > 0.5 


45 


41 


4 


24 


41 


44 


median(vpec) 


206 




149 




223 




median(vpec_i ) 


0.41 




0.28 




0.46 




Pa < 0.05 


77 


71 


7 


41 


70 


76 


Pc < 0.05 


70 


64 


5 


29 


65 


71 


Pp < 0.05 


63 


58 


1 


6 


62 


67 


Pad < 0.05 


24 


22 


3 


18 


21 


23 


Psw < 0.05 


29 


27 


3 


18 


26 


28 


Pkurtosis < 0.05 


13 


11 








13 


14 


Pskewness — 0.05 


1 


1 








1 


1 



Notes. Numbers and fractions of clusters in various populations and 
median values of the peculiar velocities of main galaxies, Vpe^, (in km/s), 
and of the normalised peculiar velocities, Vpec,r. 

Table [1] shows that more than 80% of the clusters in our 
sample have multiple components, as detected by Mclust by the 
3D analysis. About 75% of the clusters continue to show multi- 
ple components, if we use only the sky distribution of galaxies. 
More than 70% of all the clusters and more than 75% of the 
multicomponent clusters show statistically significant substruc- 
ture according to the DS test, also; these fractions are slightly 
smaller for the a and fi tests. The fractions of clusters with non- 
Gaussian velocity distributions are much smaller The reason is 
that in many clusters galaxies from different components in the 
sky coordinates have similar (nearly Gaussian) distributions of 
velocities (we will show below a few examples of that). One- 
component clusters show signatures of substructure and non- 
Gaussian velocity distributions according to the DS, a, and other 
tests as well, but as seen from Table[Tl the fractions of such clus- 
ters are smaller than those for multicomponent clusters. There 
are three clusters for which /?ad ^ 0.05, but /?sw > 0.05, and 
eight clusters with /?sw ^ 0.05, but pad ^ 0.05. iHou et al] 
(l2009h showed that the AD test is very reliable; in our study the 
SW test detected non-Gaussian velocity distributions in a larger 
number of clusters, but there are also clusters for which the AD 
test found significant non-Gaussianity in the velocity distribu- 
tions of galaxies, but the SW test did not. 



4. 1 . Peculiar velocities of cluster main galaxies 

Figure [3] shows the distributions of the peculiar velocities Vpec 
and normalised peculiar velocities Vpec,i of the main galaxy 
in clusters, separately for one-component and multicomponent 
clusters. The peculiar velocity and normalised peculiar veloc- 
ity distributions for main galaxies of one-component clusters we 
see a maximum at velocities less than 250 km/s (0.5 for the nor- 
malised peculiar velocities), followed by another maximum at 
larger velocities. Checking the distributions of galaxies in the 
clusters shows that in these clusters where the main galaxy has a 
peculiar velocity less than 250 km/s, it is located in a component 
close to the central part of the cluster If the peculiar velocity is 
higher, the main galaxy is located at the edge of the cluster, in 
multicomponent clusters typically in another component. 

Table[T]shows that in about half of the multicomponent clus- 
ters the peculiar velocity of main galaxy is larger than 250 km/s 
(0.5 for normalised peculiar velocities). There is a clear differ- 
ence between the median values of the pecuUar velocities and 
normalised peculiar velocities of the one-component and multi- 
component clusters - these velocities are much larg er in the mul- 
ticom ponent clusters. This agrees with the results of lCoziol et al.l 
(12009) who found that the median value of the normalised pe- 
culiar velocities of main galaxies in the Abell clusters is 0.32, 
which is between the values found in this study. 

0.004 I ^ , ^ , ^ 1 2.5 I ^ , ^ , ^ , ^ 1 




Vpec (km/s) Vpec.r 

Fig. 3. Distribution of the pecuUar velocities of main galaxies, 
Vpec (left panel), and the normalised peculiar velocities of main 
galaxies, Vpec.r (right panel), for clusters with one component 
(dashed line) and with multiple components (solid line). 

Figure |4] shows the 2D distance of the main galaxy from the 
cluster centre, both for the multicomponent and one-component 
clusters. It is well seen that in multicomponent clusters a large 
fraction of main galaxies are located far away from the cluster 
centre (grey dotted line). However, when looking at the compo- 
nents found by 3D normal mixture modelling, we see that the 
main galaxies of clusters are preferentially located close to the 
centre of one of the components red solid line). The distribution 
of distances from the component centre for the brightest galaxies 
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Cluster main galaxies: cluster centre 
Cluster main galaxies: component centre 
Component main galaxies: component centre 
Unimodal cluster main galaxies 
Unimodal cluster top3 galaxies 




0.5 1 
Distance from (sub)cluster centre (h'^Upc) 

Fig. 4. Distribution of the distance of the main galaxy from the 
cluster (subcluster or component) centre for various subsamples 
of galaxies. Multicomponent clusters: grey dotted line - the dis- 
tance of the main galaxy from the cluster centre; red solid line 
- the distance of the main galaxy from the (nearest) component 
centre; blue dashed line - the distance of the brightest galaxy in 
a component from the component centre. One-component clus- 
ters: dark green long-dashed line - the distance from the cluster 
centre; light green dotted-dashed line - the minimum distance of 
one of the three brightest galaxies from the cluster centre. 
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Fig. 5. Peculiar velocities of main cluster galaxies and the num- 
ber of separate components in the clusters. Large symbols mark 
clusters with a non-Gaussian velocity distribution of galaxies 
(pad < 0.05 and psw < 0.05). Empty circles denote clus- 
ters with normalised peculiar velocities of their main galaxies 
Vpec r < 0.5, filled circles denote clusters with Vpec r > 0.5. 



in the components (Fig.|4]i shows that these galaxies are also lo- 
cated preferentially close to the component centre (blue dashed 
Une). 

The distribution of the distances of the main galaxy from 
the cluster centre in one-component clusters (Fig. H) dark green 
long-dashed line) shows that more than a half of the main galax- 
ies are located near the cluster centre, but there is also a sub- 
stantial number of main galaxies, which are further away from 
the cluster centre. We calculated also the minimum distance 
from the cluster centre for three brightest galaxies in the one- 
component clusters. As seen from Fig. |4] (light green dotted- 
dashed line), one of the three brightest galaxies in clusters is 
always located close to the cluster centre. This shows that the 
central galaxy of a cluster is typically one of the most luminous 
galaxies, but not always the most luminous one. 

4.2. Comparison of the results of different tests 

In Figure|5]we compare the following properties of clusters: the 
number of components determined by the 3D version of Mclust, 
the peculiar and normalised peculiar velocities of main galax- 
ies, and the Gaussianity of the velocity distribution of galaxies 
in clusters, as estimated by the Anderson-Darling and Shapiro- 
Wilk tests. Here we see that we have clusters consisting of up 
to eight components. The peculiar velocities of the main galax- 
ies may be large both in the clusters with one component and 
with multiple components, as was seen also in Fig. [3] There are 
six clusters with five or more components, but only two of them 
have non-Gaussian velocity distributions of galaxies according 
to both tests (both are five-component clusters with Vpec ~ 280 
km/s, so their symbols coincide in Fig.|5]l. 

Figure |6]compares the results of the DS test and the 3D com- 
ponent analysis. Symbol sizes indicate the values of the peculiar 
velocities of the main galaxies in the clusters. Figure |6] shows 
that many multicomponent clusters, for which the DS test finds 
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Fig. 6. 3D tests: the Dressler-Shectman test /?-value versus the 
number of components in a cluster Symbol sizes are propor- 
tional to the peculiar velocities of main galaxies, Vpec. Small p- 
values show high probability of substructure. 



the presence of substructure, have large pecuUar velocities of 
the main galaxes. There are seven one-component clusters af- 
ter Mclust, for which the DS test suggests that these clusters 
have significant substructure (see also Table [1]). Visual analysis 
of these clusters shows that these clusters have tails in their ve- 
locity distribution, and for some of these clusters the AD and/or 
SW tests also suggest a non-Gaussian distribution of galaxy ve- 
locities. There are 22 multicomponent clusters in our sample, 
which according to the DS test have no significant substructure. 
Our analysis showed that these clusters have components in the 
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sky distribution of galaxies with similar distributions of veloci- 
ties. Thus while Mclust finds better the components in the sky 
distribution of galaxies, the DS test is more sensitive to the pos- 
sible substructure in the velocity distribution. 

Figure Q compares the results of the fi test and the 2D com- 
ponent analysis with Mclust. The maximum number of compo- 
nents detected in the sky coordinates is slightly smaller than that 
found using both the sky coordinates and the velocities of galax- 
ies. Figure |2] shows that there are multicomponent clusters with- 
out significant asymmetry {pp > 0.05) and with large peculiar 
velocities of main galaxies (30 clusters). Comparison with the 
results of the DS test shows that 20 of these clusters have sig- 
nificant substructure, according to this test. There are 12 multi- 
component clusters with significant asymmetry in the sky distri- 
bution of galaxies, according to the /3 test, but without significant 
substructure, according to the DS test, showing that the velocity 
distribution of galaxies in different cluster components is similar. 
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Fig. 7. 2D tests: the number of cluster components versus the 
asymmetry [p/j). Symbol sizes are proportional to the peculiar 
velocities of main galaxies, Vpec. 

Next we will use the principal component analysis to com- 
pare the results of all tests simultaneously. We will use logarith- 
mic scale for the peculiar velocities of main galaxies, this makes 
their range smaller. Figure [8] and Table |2] show the results of 
this analysis. In the biplot (Fig. [8]) the arrows represent the axes 
where each original variable lies, and their length is proportional 
to their importance within each PC. In calculations we use I - p 
instead of the p-value - larger values of 1 - /:> suggest a higher 
probability to have substructure, therefore in this case the arrows 
corresponding to the number of the components in 3D and 2D 
point towards the same direction as the arrows corresponding to 
the DS, a, and /3 tests. The locations of clusters in the PC space 
are also plotted. In Table |2] we give the values of the principal 
components and the standard deviations, the proportion of vari- 
ance, and the cumulative variance of the principal components. 
The values of the components show the importance of the origi- 
nal parameters in each PCi. 

In Fig. [8] the arrows corresponding to the 3D and 2D tests 
are oriented towards about the same direction, the arrows cor- 
responding to the ID tests are pointed at a direction almost at 
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Fig. 8. Principal component analysis biplot for all tests. 



the right angle to this, and form another set. This suggests that 
the results of the 3D and 2D tests are correlated, as well as the 
results of ID tests among themselves. Among the 3D and 2D 
tests the highest importance in first two principal PCA compo- 
nents have the numbers of components determined with normal 
mixture modelling, followed by the DS test. The a and /3 tests 
are of almost equal importance. Among the ID tests, the results 
of the Anderson-Darling and the Shapiro-Wilk tests have equal 
importance in determining the non-Gaussianity of the velocity 
distribution. 

Table |2] shows that seven principal components are needed 
to explain more than 90% of the variance in the parameters, thus 
there is no one dominant test which could find most of the infor- 
mation about multimodality of the galaxy distribution in clus- 
ters. 

Different types of clusters populate different regions in the 
PC1-PC2 plane. Unimodal clusters with Gaussian distribution 
of velocities are located in the upper lefthand part of the plot 
and have larger PC2 and larger negative PCI values (for ex- 
ample, the clusters 608 and 58604). Multimodal clusters popu- 
late the lower righthand area of the biplot (the clusters 34727 
and 914). Clusters with small values of the peculiar veloci- 
ties of main galaxies, with a Gaussian distribution of velocities, 
but with a large number of components populate the lefthand 
lower area of the PC1-PC2 plane (the clusters 29348 and 4744). 
Multimodal clusters with large peculiar velocities of main galax- 
ies and with non-Gaussian velocity distributions are located in 
the upper righthand area of the plane (the clusters 4122, 23374 
and others). 

4.3. Examples of multimodal and unimodal clusters 

In our sample there are eight clusters which show signatures 
of non-Gaussianity in their galaxy distribution by all the tests. 
They have multiple components, the 3D and 2D tests confirm 
the presence of substructures and of asymmetry, and their ve- 
locity distributions differ from Gaussian at a very high statisti- 
cal significance: the clusters 880, 4122, 28387, 34726, 58305, 
67297, 68625, and 73088. Many clusters with multiple compo- 
nents are missing from this list, for example the cluster 914, a 
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Table 2. Results of the principal component analysis for the tests. 





PCI 


PC2 


PC3 


PC4 


PC5 


PC6 


PC7 


PC8 


PC9 


PC 10 


PCll 


^corap3D 


0.218 


-0.481 


-0.098 


-0.249 


-0.016 


0.363 


-0.060 


0.680 


-0.215 


0.040 


0.025 


^comp2D 


0.158 


-0.485 


-0.199 


-0.309 


0.011 


0.324 


-0.154 


-0.652 


0.217 


-0.007 


-0.020 




0.205 


-0.341 


0.112 


0.553 


-0.070 


0.127 


0.698 


-0.103 


-0.028 


-0.037 


-0.003 


Pa 


0.136 


-0.342 


0.158 


0.589 


0.175 


-0.210 


-0.642 


-0.001 


-0.031 


0.047 


0.031 


PP 


0.179 


-0.332 


-0.115 


-0.313 


0.218 


-0.798 


0.214 


0.012 


-0.065 


-0.086 


-0.051 


log(Vpec) 


0.301 


0.222 


-0.544 


0.184 


0.101 


0.072 


-0.071 


0.051 


-0.001 


-0.706 


-0.076 




0.242 


0.185 


-0.624 


0.158 


0.062 


-0.034 


0.031 


-0.044 


-0.145 


0.678 


0.048 


Pad 


0.499 


0.169 


0.195 


-0.085 


0.018 


-0.026 


0.013 


0.084 


0.438 


0.022 


0.689 




0.494 


0.170 


0.239 


-0.053 


-0.035 


0.011 


-0.030 


0.096 


0.352 


0.149 


-0.712 


Pkiiit 


0.297 


0.035 


0.071 


-0.051 


-0.812 


-0.161 


-0.135 


-0.151 


-0.410 


-0.065 


0.048 


Pskew 


0.320 


0.211 


0.334 


-0.136 


0.490 


0.177 


0.019 


-0.234 


-0.627 


-0.022 


0.026 


Importance of the components 




PCI 


PC2 


PC3 


PC4 


PC5 


PC6 


PC7 


PC8 


PC9 


PC 10 


PCll 


St. deviation 


1.783 


1.392 


1.283 


1.090 


0.969 


0.863 


0.691 


0.649 


0.524 


0.370 


0.227 


Prop, of Variance 


0.289 


0.176 


0.150 


0.108 


0.085 


0.067 


0.043 


0.038 


0.025 


0.012 


0.004 


Cumulative Prop. 


0.289 


0.465 


0.615 


0.723 


0.808 


0.876 


0.919 


0.957 


0.982 


0.995 


1.000 



Notes. Last three rows present the standard deviation, proportion of variance, and cumulative variance of principal components. 
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Fig. 9. Two of the most non-Gaussian clusters in our sample. Upper row: the cluster 34726 (Abell cluster A2028). Lower row: the 
cluster 914 (A1750). From left to right: the DS test bubble plot (here symbol sizes are proportional to the e^), and R.A. vs. Dec, 
and R.A. vs. velocity (in Kfikml sec) plots; the symbols show different components as found by Mclust. The star marks the location 
of the main galaxy. 



well-know n binary merging A bell cluster A1750 in the Sloan 
Great Wall (lEinasto et alj20ll|201]ilS, and references therein), 
because in many multicomponent clusters galaxies from differ- 
ent components have similar distributions of velocities. 

Figure |9] shows the distribution of galaxies in a cluster with 
four components, the cluster 34726, and in the cluster 914 with 
five components. 

The cluster 34726 is one of the two clusters in our sam- 
ple which can be identified with the Abell cluster A2028 in 
sky projection, another one is the clust er 34727. A2028 is a 
merging X-ray cluster with substructure (IFlin & Krvwultll2006t 



iGastaldello et alJbOlOl) . which redshift close to the redshift of 
our cluster 34727. This is one example of projections in the 
Abell cluster catalogue. Figure |9] (upper row) shows that the sky 
distribution of galaxies in the cluster 34726 has two main com- 
ponents. In one of them the concentration of galaxies is espe- 
cially high. This component is well seen in the distribution of 
velocities where very large velocity dispersion suggest a finger- 
of-god effect (right panel). Here also according to the DS test 
the probability of substructure is high. The a and fi tests confirm 
the presence of substructure and asymmetry in the distribution 
of galaxies, the AD and SW tests show non-Gaussianity of the 
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Fig. 10. Two of the most Gaussian clusters in our sample. Upper row: the cluster 58604. Lower row: the cluster 60539. From left to 
right: the DS test bubble plot (symbol sizes are proportional to the e*), and R.A. vs. Dec, and R.A. vs velocity (in lO^km/s) plots; 
the symbols show different components as found by Mclust. The star marks the location of the main galaxy. 



distribution of velocities of galaxies in this cluster The pecuUar 
velocity of the main galaxy is small. 

The cluster 914 (Fig. |9] lower row) that corresponds to the 
Abell c luster A 1750 is a well-known merging bi nary X-ray 
cluster (iDonnellv et al. 2001; Belsole et al. 2004; Burgett et al.1 
I2004t iHwang & Le^009; .Einasto et al.,,2010 ). In this cluster 



we see even five components. The DS test shows that the prob- 
ability of substructure is the largest in the component where the 
main galaxy is located (Fig. |9] lower left panel). The peculiar 
velocity of the main galaxy is large. Other 3D and 2D tests con- 
firm the presence of substructure and an asymmetrical galaxy 
distribution in A 1750. Galaxies from different components in 
this cluster have similar velocity distributions, so according to 
the AD and SW tests their dis tribution is Gaussia n (for details 
and references we refer also to lEinasto et al.ll20ldh . 

Next we list the clusters for which all tests confirm uni- 
modality and the Gaussian velocity distribution of galaxies: 608, 
5217, 25078, 39914, 50129, 58604, and 67116. In Fig.[l0]we 
plot the distribution of galaxies in two of them - in the cluster 
58604, and in the cluster 60539 (A1516, the richest cluster in 
the superclus ter SCI 1 1 1 in the Sloan Great Wall, described in 
lEinasto et al.l ^OlO). For the cluster 58604 all our tests confirm 
that this is an unimodal cluster. The peculiar velocity of the main 
galaxy is rather large for an unimodal cluster, about 200 km/s. 
The cluster A1516 has the most regular ellipsoidal sky distribu- 
tion of galaxies among the clusters in our sample, with a small 
peculiar velocity of the main galaxy. However, the distribution 
of velocities in the central part of the cluster shows an excess of 
negative (in respect of the cluster center) velocities, and accord- 
ing to the DS and a tests this cluster may have substructure (see 
also an analysis of this cluster in Einasto et al. 201 0). This may 
be a s ignature of a line-of-sight merger (see also [Pinknev et al.l 
II996h . 



In Tables lA.ll - IA.4I we present the data about the most 
unimodal and most multimodal clusters, as well as data about 
the clusters with at least five components, and about those one- 
component clusters for which the DS test found significant sub- 
structure. 



4.4. Non-Gaussianity parameters and physical properties of 
clusters 

The principal component analysis can help us to study the re- 
lation between the physical parameters of clusters and the pres- 
ence of substructure. In these calculations we include the num- 
ber of components as determined with the 3D normal mixture 
modelling, the /?-value for the DS test, p/^, the peculiar velocity 
of the main galaxy in a cluster, Vpec, the number of galaxies in a 
cluster, A^gai, the total luminosity of a cluster, Ltot, the virial ra- 
dius of a cluster, rvh , and the rms velocity of the cluster galaxies. 
We again use 1 - pA, and use logarithms of the peculiar veloc- 
ities of main galaxies and of the physical parameters. Figure [TTI 
and Table|3]show the results of this analysis. 

In Fig.[TT]the arrows corresponding to the tests for substruc- 
ture and the arrows corresponding to the physical parameters of 
clusters do not point to the same direction - the correlations be- 
tween the parameters are not strong, when we consider all pa- 
rameters simultaneously. The coefficients of the first principal 
component in Table [3] show that the physical parameters are of 
larger importance than the substructure indicators. The coeffi- 
cient corresponding to the luminosity of a cluster is the largest, 
but the differences between the values of the coefficients are not 
large, so there is no single dominant parameter Among the co- 
efficients of the second principal component the number of com- 
ponents has the largest (absolute) coefficient. Among the coef- 
ficients of the third principal component the peculiar velocities 
of the main galaxy, and the rms velocity of galaxies in a clus- 
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Table 3. Results of the principal component analysis for the test results and for the physical parameters of clusters, combined. 





PCI 


PC2 


PC3 


PC4 


PC5 


PC6 


PC7 


^corap3D 


0.218 


-0.596 


0.407 


-0.115 


0.600 


0.232 


-0.050 


log(Vpec) 


0.185 


0.305 


0.415 


-0.811 


-0.147 


-0.133 


0.039 


PA 


0.284 


-0.171 


0.602 


0.427 


-0.571 


-0.051 


0.118 


log(L,oi) 


0.537 


0.116 


-0.287 


-0.010 


0.035 


0.351 


0.699 


log(iVg,i) 


0.509 


-0.004 


-0.153 


0.120 


0.242 


-0.799 


-0.063 


log{rms velocity) 


0.355 


0.603 


0.162 


0.268 


0.236 


0.343 


-0.487 


log(rvir) 


0.401 


-0.379 


-0.406 


-0.242 


-0.418 


0.211 


-0.500 


Importance of components 




PCI 


PC2 


PC3 


PC4 


PC5 


PC6 


PC7 


St. deviation 


1.627 


1.115 


1.014 


0.950 


0.795 


0.609 


0.418 


Prop, of Variance 


0.378 


0.177 


0.147 


0.129 


0.090 


0.053 


0.025 


Cumulative Prop. 


0.378 


0.556 


0.703 


0.832 


0.922 


0.975 


1.000 



ter are of equal importance. In Fig. [TT] the arrow corresponding 
to the peculiar velocity of the main galaxy in a cluster points to 
the same direction as the arrow for the rms velocity of a cluster, 
showing that these velocities are larger in clusters with larger 
rms velocities. These clusters are also richer. 
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Fig. 11. Principal component analysis with A^comp3D, Vpec, and 
Pa, and the physical parameters of clusters (L,o,, Ngai, fvir, and 
rms velocity). 

TableOshows that the first principal component accounts for 
almost 38% of the variance of parameters. Five principal com- 
ponents are needed to explain more than 90% of the variance of 
the parameters - clusters are complicated objects, the properties 
of which cannot be explained with a small number of parameters 
(in contrast to supercluster s which can be descr ibed with a small 
number of parameters, see Einasto et ^ 1201 lah . 

In the PC1-PC2 plane one-component clusters with a 
Gaussian distribution of velocities are located at the upper left- 
hand part of the plot and have larger PC2 and smaller (larger 
negative) PCI values (for example, the clusters 608 and 28508). 
Multicomponent clusters of high luminosity populate the lower 
and middle righthand area of the biplot (the clusters 34727 and 
914). Less luminous multicomponent clusters populate the left- 
hand lower area of the PCI- PC2 plane (the clusters 1 1474 and 



11015). One-component clusters with large rms velocities of 
galaxies populate the upper righthand area of the plane (the clus- 
ters 16350, 17210, and others). 

We checked for pairwise correlations between the physical 
parameters of clusters and the substructure characteristics and 
found that correlations between the number of galaxies in clus- 
ters, the total luminosity of clusters, the virial radius of clusters, 
and the numbers of components and the presence of substructure 
according to the 3D and 2D tests are statistically highly signif- 
icant - richer, larger and more luminous clusters have a larger 
amount of substructure. 



5. Discussion and conclusions 

5.1. Comparison witli otiier studies 

An increasing number of studies have shown recently the pres- 
ence of substructure in groups and clusters of galaxies deter- 
mined using optical or X-ray data (see references in Sect.[T]). We 
detected multiple components in more than 80% of the clusters 
using multidimensional normal mixture modelling. Comparison 
with other studies in Sect.|4]showed that although this is a higher 
fraction of multimodal clusters than found in other studies, com- 
parison of the fractions of clusters with substructure or with non- 
Gaussian velocity distributions determined with similar methods 
gives the results in a good agreement with others. This supports 
our choice of the parameters for the FoF algorithm for group 
definition and suggests that in our catalogue groups with sub- 
structure are real groups and not complexes of small groups, ar- 
tificially linked together by an unreasonable choise of linking 
parameters. 

The fra ctions of clusters with substructure found in X-ray 
clusters by iRamella et alJ (l2007h is almost the same that we 
found with the 2D normal mixture modelling. In these studies 
there are eight common clusters; i n seven of the m substructure 
have been detected in both studies. iRamella etal.l (12007) found 
two components in A602 (cluster 43545), this is unimodal clus- 
ter according to our calculations. In the sky distribution this clus- 
ter has a very elongated shape; this may be visible in X-rays as 
two components. 

Comparison with lFlin & KrvwuE (12006 ) shows that they de- 
tected substructures in all 14 clusters common for both studies; 
we found multiple components in 12 clusters. A1809 (our clus- 
ter 671 16) is an unimodal cluster according to Mclust, although 
the DS and a tests detected substructure at the 90% signifi- 
cance level. We detecte d no substructure in A1650 (25078) while 
lFlin&Krvwuia(l2006h mention that there is one group of galax- 
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ies in the field of the galaxies in this cluster. Ther efore, an overall 
agreement between the results of our study and lFlin & Krvwuli 
(2006) is very good, considering that the sample selection and 
the methods to analyse modality of clusters are different. 

The cluster A1750 (our cluster 914) has been studied for 
substructure by other authors, too. We mentioned above that 
the components found in this cluster by us coincide well with 
those found in other studies (see Einast o et a l. 2010, for details). 
We foun d an especially good agre ement between our re sults and 
those bv lBurgett et"alT(l2004 and lHwang & Leel(l2009h . 

We hav e three common clus ters with those searched for 
rotation by Hwang & Led (l2007h . We found that all of them 
(A1035, A1139, and A2169; these are our clusters 20159, 3714, 
and 12508) are multi modal clusters which agrees with the re- 
sults of the DS test bvl Hwang & Led (l2007l) . They also suggest 
that the clusters A1139 and A2169 may be rotating or merg- 
ing. We have also fi ve common clusters with those studied by 
iBurgett et all ( |2004 using data from the 2dFGRS. All of them 
(Al 139, A1238, A1620, A1663, and A1750, these are our clus- 
ters 3714, 4744, 24554, 24829, and 914) are multimodal clus- 
ters, as found in b oth studies. Eight of te n common clusters with 
those studied by lOegerle & Hilll (1200 ll) are multimodal in our 
study, for one unimodal cluster (A2089, our cluster 39914) the 
Pa = 0.05 1 tha t show s the presence of substructure, as found in 
lOegerle &Hiil(l200Th . but not at a very high significance level. 

We showed that the peculiar velocities of main galaxies in a 
large fraction of clusters are large, and the most luminous galaxy 
(cluster main galaxy) is often not located close to the cluster cen- 
tre. This means that the most luminous galaxies in clusters are 
not the central galaxies. Actually, there is no stringent reason 
for that. In the merger scenario, the clusters form through merg- 
ing of smaller groups/clusters. The main galaxy in a new cluster 
is on e of the main galaxie s of the merged clusters. As pointed 
out in lTempel et alJ (120091) . the second luminous galaxies in rich 
clusters have been main galaxies before the latest merger event 
between the clusters. In this study, we analyse the richest clus- 
ters in the SDSS sample, and such rich clusters form through 
merging of other clusters. The presence of substructure in these 
clusters supports this idea. For the main galaxies this means that 
the main galaxies in rich clusters have not yet found their place 
in the cluster centre, they are still located close to the centres 
of their parent clusters (subclusters of the new cluster). Large 
distances from the cluster centres and large peculiar velocities 
of the brightest galaxies in clusters indicate that clusters are not 
yet fully virialized after the last merge r Similar results were ob- 
tained by lEinasto et alJ (l2010t l201ld) in the study of groups in 
the Sloan Great Wall. 

Cosmological simulations show the merging and growth of 
dark matter haloes (Richstone e t al. 1992; Mo & White 2002; 
iMcIntosh et aP [2008: Fakhouri etalJ 120101: iPower et alJ l201 1 . 
and references therein). The late time formation of the main 
haloes and the number of re cent major mergers can cause the late 
time subgrou ping of halo es dSmith & Tavlod2008t lEinasto et al.l 
I20r0- Power et al.ll201ll) . 

The principal component analysis of the substructure indi- 
cators suggests that the results of the 3D and 2D tests are corre- 
lated, as well as the results of the ID tests. Among the 3D and 2D 
tests the highest importance have the 3D normal mixture mod- 
elling and the DS test. With the 3D tests we detected the largest 
fraction of clusters with substructure, followed by the 2D tests. 
The PCA showed that there is no single dominant test which 
could find most of the information about multimodality of the 
galaxy distribution in clusters. Also Pin knev et al.l (Il996h con- 
cluded that the higher the dimensionality of the test, the higher 



is the probability to detect substructure. However, the tests are 
sensitive to the different aspects of substructure and to non- 
Gaussianity, thus it is preferable to use several tests of different 
dimensionality. 

Richer, larger and more luminous clusters have a larger 
amount of substructure. The principal component analysis us- 
ing both the substructure indicators and the physical parameters 
of clusters showed that five principal components are needed to 
explain more than 90% of variance in data - galaxy clusters are 
complicated objects, the properties of which cannot be explained 
with a small number of parameters, as was shown recently also 
bv lJeeson-Daniel etali j201lb and ISkibba & MacciSi (l201lb by 
the PCA analysis of dark matter haloes. 

Cosmological simulations show that the fraction of sub - 
structure is higher in h igh-r edshift ha l oes ( Giocoli et alJ 1201 Ol) . 
ISwinbank etall (l2007l) and iGal et all (l2008h found evidence of 
substructure in clusters of high-redshift superclusters at z » 0.9. 
More data about high-redshift clusters is needed to make a sta- 
tistical comparison. 

Current theories of the formation and evolution of galaxies, 
and galaxy groups and clusters tell that galaxi es and their sys- 
tems form in vi rialised dark matter haloes (Whi te & Reesl[T978l : 
White & Frenk 1991), which grow hierarchically by merging of 
smaller mass haloes (see, e.g., L oeb 2008 ). The halo model of 
group and cluster properties and clustering has been success- 
ful in explaining several properties of groups and clusters, their 
galaxy population, and clusterin g (Berlind & Weinberg 20()3; 
iSer hnd et alJ2006Hvan den Bosch et al. 2007; Skibba 2009, and 
references therein). The halo model h as been also used i n 
studies of the structure of superclusters (lEinasto et alJl20lT^ . 
In virialised clusters galaxies follow the gravitational poten- 
tial and the main galaxies of clusters (their brightest galax- 
ies) lie at the centers of clus t ers and have small peculiar ve- 
lociti es jOstriker & Tremaind 11 9751: iMerrittI 119841: iMalumuthI 
[1991 [lerlind & Weinberg 2002^ \Yi 



"ang et al 



2005). The most 



important halo parameter is the halo mass, although several 
studies have shown that other parameters are also impor- 
tant in shaping the properties of dark ma tter haloes, like the 
formatio n time, concentration and others ( Jeeson-Daniel et aP 
2011; S kibba & Maccioll201 lb ICroft et al.l[2oT l], and references 
therein). 

We showed that a significant fraction of the main galaxies 
in clusters have large peculiar velocities, and more than 80% of 
clusters in our sample have substructure and/or non-Gaussian 
velocity distributions. Knebe & Miiller (2000) applied the DS 
statistics to search for substructure of clusters in numerical sim- 
ulations. They studied clusters from a set of cosmological mod- 
els and found that in all models, the DS test detected significant 
substructure in 30-40 % of clusters, which is about half of what 
we found in this study. The high fraction of clusters with sub- 
structure and large peculiar velocities of main galaxies show that 
clusters are not in dynamical equilibrium. We also found that 
galaxy clusters are complicated objects, the properties of which 
cannot be explained with a small number of parameters. These 
differences between observations and simulations do not fit well 
into the halo model framework and have yet to be explained. 



5.2. Conclusions 

We searched for substructure and non-Gaussian velocity distri- 
butions in rich clusters drawn from the SDSS DR8. We present 
lists of unimodal and multimodal clusters. Our conclusions are 
as follows. 
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1) We showed, using a number of tests, that more than 80% of 
rich clusters have substructure and/or non-Gaussian velocity 
distributions of galaxies. 

2) The peculiar velocities of the main galaxies in clusters, and 
their distances from the cluster centre are large, especially in 
clusters with multiple components. In multicomponent clus- 
ters the brightest galaxies are typically located near the cen- 
tre of one of the components. 

3) The largest number of clusters with substructure was de- 
tected by multidimensional normal mixture modelling, fol- 
lowed by the Dressler-Shectman test. 

4) Principal component analysis shows that there is no one 
dominant test which could find most of the information about 
the multimodality of galaxy distribution in clusters. Different 
tests are sensitive to different signatures of multimodality, 
therefore it is important to use several tests to search for sub- 
structure in clusters. 

5) Richer, larger, and more luminous clusters with larger rms 
velocities have larger amount of substructure and larger (nor- 
malised) peculiar velocities of the main galaxies. 

6) Principal component analysis using both substructure indi- 
cators and the physical parameters of clusters showed that 
galaxy clusters are complicated objects, the properties of 
which cannot be explained with a small number of param- 
eters. 

7) Our results show that simple halo model do not explain all 
the properties of observed clusters. The halo model assumes 
that haloes (clusters) are virialised, while we found that they 
are not. Also, the fraction of observed clusters with substruc- 
ture is larger than that found in simulations. 

The presence of substructure, large distances of main galax- 
ies from the cluster centre, and their large peculiar velocities is a 
sign of mergers and/or infall, and suggests that most clusters in 
our sample are not yet in dynanoical equilibrium. The high fre- 
quency of such clusters tells that mergers between groups and 
clusters are common - galaxy groups continue to grow and are 
still assembling. Our unimodal clusters are examples of clusters 
which are probably already in dynamical equilibrium and there- 
fore can be used in the studies of cluster characteristics Uke mass 
estimation, analysis of concentration and others reliably. 

To understand better the properties of galaxy clusters and 
hence the formation and evolution of the structures in the 
Universe we plan to analyse in detail galaxy populations in dif- 
ferent cluster components, and the connection of the properties 
of clusters and the environment where they reside. 
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Table A.l. Data on unimodal clusters 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


ID 




R.A. 


Dec. 


Dist. 








Abell ID 






[deg] 


[deg] 


[h-^ Mpc] 


[lo'^r^Lo] 




[r' Mpc] 




608 


60 


245.23 


29.83 


284.84 


132.16 


532.83 


0.66 


A 2175 


748 


79 


159.85 


5.22 


206.78 


93.59 


748.24 


0.43 


A 1066 


5217 


89 


180.06 


56.22 


193.33 


94.88 


577.08 


0.61 


A 1436 


25078 


51 


194.69 


-1.70 


249.24 


87.64 


498.42 


0.61 


A 1650 


28272 


51 


248.33 


11.80 


153.80 


40.89 


355.14 


0.50 


— 


32663 


51 


225.59 


21.36 


186.20 


61.37 


435.82 


0.73 


— 


39914 


63 


233.12 


28.02 


218.38 


69.77 


446.50 


0.65 


A 2089 


43336 


68 


207.98 


46.33 


186.81 


73.08 


471.40 


0.60 




43545 


51 


118.25 


29.41 


179.41 


42.75 


577.85 


0.51 


A 602 


50129 


52 


205.43 


26.39 


225.02 


61.83 


449.92 


0.51 


A 1775 


58604 


58 


170.39 


2.86 


150.50 


39.84 


528.86 


0.42 




60539 


107 


184.40 


3.65 


230.10 


136.74 


830.47 


0.55 


A 1516 


63361 


72 


205.59 


2.26 


230.57 


103.06 


646.10 


0.64 


A 1773 


67116 


80 


208.29 


5.17 


236.99 


114.14 


651.61 


0.44 


A 1809 



Notes. Columns are as follows: 1: ID of the cluster; 2: the number of galaxies in the cluster, Afgai; 3^: cluster right ascension and declination; 5: 
cluster comoving distance; 6: cluster total luminosity; 7: rms velocity of galaxies in the cluster; 8: cluster virial radius; 9: Abell ID of the cluster. 



Table A.2. Results of the tests. Unimodal clusters 



(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


ID 








A^cSD 




uriiD 


Pa 


Pa 


Nc2D 


urilD 




Pad 


Psw 


Pkun 


/-'skew 


608 


60 


-228.43 


-0.429 




1 


•10-^ 


0.884 


0.080 




1 ■ 10-' 


0.299 


0.735 


0.860 


0.611 


0.511 


748 


79 


613.95 


0.821 




1 


•10-3 


0.007 


0.007 




1 • 10-3 


0.011 


0.308 


0.076 


0.154 


0.564 


5217 


89 


-19.11 


-0.033 




1 


• 10-3 


0.216 


0.110 




1 • 10-3 


0.367 


0.240 


0.309 


0.136 


0.658 


25078 


51 


-120.45 


-0.242 




1 


•10-3 


0.270 


0.378 




1 • 10-3 


0.829 


0.258 


0.513 


0.753 


0.452 


28272 


51 


464.63 


1.308 




1 


•10-3 


0.032 


0.213 




1 • 10-3 


0.645 


0.244 


0.307 


0.568 


0.298 


32663 


51 


96.03 


0.220 




1 


•10-3 


0.008 


0.147 




1 • 10-3 


0.564 


0.060 


0.041 


0.174 


0.118 


39914 


63 


148.99 


0.334 




1 


•10-3 


0.051 


0.106 




1 • 10-3 


0.804 


0.904 


0.863 


0.566 


0.759 


43336 


68 


-97.43 


-0.207 




1 


■ 10-3 


0.003 


1 ■ 10-3 




1 • 10-3 


0.335 


0.980 


0.941 


0.756 


0.782 


43545 


51 


81.55 


0.141 




1 


• 10-3 


0.004 


0.011 




1 • 10-3 


0.055 


0.047 


0.200 


0.817 


0.369 


50129 


52 


-129.26 


-0.287 




1 


•10-3 


0.326 


0.737 




1 • 10-3 


0.357 


0.446 


0.421 


0.181 


0.770 


58604 


58 


-201.00 


-0.380 




1 


•10-3 


0.504 


0.232 




1 • 10-3 


0.222 


0.522 


0.677 


0.363 


0.847 


60539 


107 


-24.93 


-0.030 




1 


•10-3 


0.021 


0.002 




1 • 10-3 


0.563 


0.776 


0.645 


0.882 


0.479 


63361 


72 


-175.06 


-0.271 




1 


•10-3 


0.004 


0.179 




1 • 10-3 


0.278 


0.240 


0.455 


0.347 


0.820 


67116 


80 


-183.37 


-0.281 




1 


•10-3 


0.094 


0.089 


2 


0.149 


0.342 


0.929 


0.805 


0.930 


0.862 



Notes. Columns are as follows: 1: ID of a cluster; 2: the number of galaxies in a cluster, Wgai; 3: the peculiar velocity of the main galaxy (fa/i.s--'), 
4: the normalised peculiar velocity of the main galaxy, v^^lcr^; 5: the number of components in 3D in a cluster, A'^omp; 6: the median uncertainty 
of classification (hereafter 1 • 10-3 denotes the value < 1 • 10-3); 7: the p-value for the DS test; 8: the p-value for the a test; 9: the number of 
components in 2D in a cluster, A'comp; 10: the median uncertainty of classification in 2D; 11: the p-value for the jS test; 12: the p-value for the AD 
test; 13: the p- value of the SW test; 14: the p- value for the kurtosis test; 14: the p- value for the skewness test. 
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Table A.3. Data on multimodal clusters 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


ID 




R.A. 


Dec. 


Dist. 




cr 


'"vir 


Abell ID 






[deg] 


[deg] 


Mpc] 


[lo'^r^Lo] 




[h-^ Mpc] 




880 


57 


211.40 


6.27 


249.48 


101.82 


411.37 


0.84 


— 


914 


119 


202.65 


-1.94 


256.38 


227.45 


657.41 


0.83 


A 1750 


4122 


88 


173.06 


56.09 


159.40 


68.62 


963.42 


0.49 


A 1291 


11015 


52 


212.51 


55.07 


124.69 


45.80 


303.10 


0.46 


— 


11474 


51 


218.27 


52.90 


134.86 


36.44 


306.38 


0.49 


— 


28387 


88 


169.20 


54.53 


208.84 


121.26 


481.54 


0.66 


— 


28986 


66 


233.43 


31.07 


197.06 


73.09 


398.63 


0.69 


A 2092 


34726 


145 


228.09 


7.44 


135.44 


121.20 


506.15 


0.74 


A 2028,A 2033,A 2040 


34727 


256 


227.74 


5.78 


234.02 


351.58 


825.94 


1.25 


A 2028,A 2029,A 2033,A 2040 


58305 


167 


223.28 


16.76 


135.12 


120.04 


401.84 


0.61 


A 1983 


62138 


124 


223.64 


18.67 


175.90 


127.48 


456.42 


0.78 


A 1991 


67297 


95 


127.24 


30.48 


150.20 


86.98 


770.39 


0.45 


A 671 


68625 


92 


231.04 


29.90 


336.50 


301.28 


874.84 


0.79 


A 2069 


73088 


141 


215.50 


48.40 


213.15 


184.99 


631.93 


0.71 


A 1904 



Notes. Columns are as in Table lA.l] 



Table A.4. Results of the tests. Multimodal clusters 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


ID 






Vpecr 


A'cSD 


uniD 


Pa 


Pa 


A'c2D 


un2D 


Pa 


Pad 


Psw 


Pkmt 


Pskew 


880 


57 


400.67 


0.974 


3 


1 • 10--' 


0.002 


0.002 


2 


0.001 


0.043 


0.002 


0.006 


0.135 


0.078 


914 


119 


-704.39 


-1.071 


5 


0.010 


4 • 10-^ 


0.001 


5 


0.001 


1 • 10-3 


0.140 


0.189 


0.297 


0.475 


4122 


88 


-1091.85 


-1.133 


3 


0.002 


4 • lO-'^ 


1 • w-^ 


2 


1 • 10-3 


1 • 10-3 


1 • 10-3 


1 ■ 10-3 


1 • 10-3 


0.129 


11015 


52 


140.80 


0.465 


8 


1 ■ 10-3 


0.037 


0.232 


3 


1 • 10-3 


0.003 


0.497 


0.627 


0.307 


0.909 


11474 


51 


-267.37 


-0.873 


5 


1 • IQ-^ 


0.001 


0.117 


4 


0.002 


1 • 10-3 


1 • 10-3 


0.004 


0.072 


0.140 


28387 


88 


-167.56 


-0.348 


4 


0.002 


4 • lO-'^ 


0.002 


3 


0.001 


1 • 10-3 


0.001 


0.017 


0.217 


0.469 


28986 


66 


-205.98 


-0.517 


5 


0.003 


0.002 


0.002 


4 


0.001 


0.003 


0.255 


0.344 


0.397 


0.929 


34726 


145 


-56.63 


-0.112 


4 


0.015 


0.001 


0.034 


2 


1 • 10-3 


1 • 10-3 


0.005 


0.016 


0.262 


0.242 


34727 


256 


-290.71 


-0.352 


5 


0.008 


4 • lO-'' 


0.008 


6 


0.017 


1 • 10-3 


0.042 


0.071 


0.120 


0.337 


58305 


167 


115.32 


0.287 


3 


0.003 


4 • 10-5 


0.017 


4 


0.009 


1 • 10-3 


1 • 10-3 


1 • 10-3 


0.010 


0.007 


62138 


124 


118.96 


0.261 


6 


0.017 


4 • 10-5 


0.030 


4 


0.046 


0.392 


0.238 


0.189 


0.912 


0.754 


67297 


95 


123.06 


0.160 


2 


0.022 


0.001 


0.003 


2 


0.106 


0.003 


0.002 


0.001 


0.356 


0.116 


68625 


92 


825.69 


0.944 


3 


0.001 


4 • 10-5 


1 • lo-"* 


2 


1 • 10-3 


0.040 


0.003 


0.015 


0.017 


0.427 


73088 


141 


-224.03 


-0.355 


3 


0.014 


4 • 10-5 


1 • 10-3 


4 


0.044 


1 • 10-3 


0.005 


0.005 


0.906 


0.113 



Notes. Columns are as in Table lA.2l The p-value 4 • 10 5 denotes the value p < 4 • 10 
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Table 1. Data on clusters 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8j 


(9) 


ID 




R.A. 


Dec. 


Dist. 






'"vir 


Abell ID 






[deg] 


[deg] 


[A"' Mpc] 


[lo'^r^L©] 


ykms \ 


Mpc] 




18 


87 


163.28 


55.04 


216.03 


110.83 


513.92 


0.73 


— 


323 


67 


241.45 


33.34 


178.18 


73.10 


276.83 


0.67 


— 


608 


60 


245.23 


29.83 


284.84 


132.16 


532.83 


0.66 


2175 


748 


79 


159.85 


5.22 


206.78 


93.59 


lA'&.lA 


0.43 


1066 


793 


122 


234.96 


21.73 


123.83 


63.65 


515.27 


0.56 


2107 


880 


57 


211.40 


6.27 


249.48 


101.82 


411.37 


0.84 


— 


914 


119 


202.65 


-1.94 


256.38 


227.45 


657.41 


0.83 


1750 


1469 


56 


151.94 


0.53 


287.04 


99.79 


418.98 


0.69 


933 


1573 


57 


153.48 


-0.89 


136.69 


35.71 


744.00 


0.25 


— 


1944 


60 


239.76 


18.06 


137.41 


47.99 


440.72 


0.46 


— 


2067 


62 


155.48 


38.53 


163.84 


66.62 


574.01 


0.44 


— 


3714 


82 


164.58 


1.56 


121.07 


54.71 


344.93 


0.60 


1139 


4122 


88 


173.06 


56.09 


159.40 


68.62 


963.42 


0.49 


1291 


4713 


80 


176.84 


55.69 


155.28 


80.31 


637.82 


0.43 


1377 


4744 


71 


170.78 


1.05 


221.36 


88.70 


492.96 


0.62 


1238 


4992 


68 


228.80 


4.37 


288.03 


157.59 


619.94 


0.65 


2048 


5217 


89 


180.06 


56.22 


193.33 


94.88 


577.08 


0.61 


1436 


7102 


54 


222.21 


11.27 


156.93 


33.32 


465.99 


0.46 


— 


7932 


50 


134.55 


38.49 


275.41 


99.94 


413.99 


0.55 


724 


9029 


78 


146.58 


43.18 


216.96 


99.94 


312.14 


0.72 


865 


9350 


89 


230.41 


7.71 


134.14 


55.05 


480.48 


0.65 


2055, 2063 


9985 


58 


152.53 


54.45 


138.94 


41.58 


384.99 


0.43 


— 


10438 


65 


158.30 


56.83 


136.38 


45.02 


395.18 


0.46 


— 


11015 


52 


212.51 


55.07 


124.69 


45.80 


303.10 


0.46 


— 


11474 


51 


218.27 


52.90 


134.86 


36.44 


306.38 


0.49 


— 


11683 


54 


183.64 


59.92 


179.36 


52.25 


342.89 


0.46 


1507 


12508 


93 


243.68 


49.32 


171.06 


90.76 


385.94 


0.77 


2169 


12540 


103 


204.00 


59.18 


209.84 


107.55 


764.36 


0.52 


1767 


13216 


57 


118.98 


45.64 


155.75 


48.58 


400.40 


0.54 


— 


13347 


50 


156.40 


47.59 


186.67 


64.31 


479.41 


0.48 


1003 


13408 


58 


240.25 


54.01 


192.66 


62.67 


282.97 


0.73 


2149 


16094 


71 


197.79 


39.26 


216.08 


113.02 


735.05 


0.50 


1691 


16309 


69 


255.66 


34.03 


290.56 


151.63 


878.96 


0.46 


2244 


16350 


65 


255.67 


33.50 


260.44 


135.10 


859.11 


0.36 


2245 


17210 


72 


257.43 


34.48 


246.34 


120.38 


872.61 


0.47 


2249 


18029 


53 


205.69 


29.94 


128.59 


37.05 


597.71 


0.41 


1781 


18048 


78 


173.28 


14.42 


242.21 


121.96 


596.51 


0.61 


— 


20159 


52 


158.02 


40.16 


203.31 


59.80 


517.98 


0.36 


1026, 1035 


20419 


58 


202.31 


37.53 


169.15 


44.07 


424.97 


0.59 


1749 


20514 


56 


168.82 


25.83 


145.55 


33.03 


317.82 


0.42 


— 


21573 


50 


176.47 


15.50 


205.25 


60.55 


364.05 


0.49 


— 


22572 


77 


167.10 


44.08 


176.56 


75.43 


533.27 


0.43 


1169 


23374 


114 


207.22 


26.68 


186.81 


100.95 


662.05 


0.72 


1795, 1818 


23524 


50 


207.65 


29.31 


228.53 


63.09 


304.50 


0.67 


— 


24554 


50 


192.37 


-1.63 


250.75 


83.48 


658.24 


0.58 


1620 


24604 


50 


209.82 


28.01 


222.60 


59.36 


857.98 


0.50 


1831 


24829 


77 


195.63 


-2.57 


246.69 


126.42 


534.45 


0.66 


1663 


25078 


51 


194.69 


-1.70 


249.24 


87.64 


498.42 


0.61 


1650 


28272 


51 


248.33 


11.80 


153.80 


40.89 


355.14 


0.50 


— 


28387 


88 


169.20 


54.53 


208.84 


121.26 


481.54 


0.66 


— 


28508 


58 


172.42 


54.10 


205.02 


74.36 


494.72 


0.40 


1270 


28986 


66 


233.43 


31.07 


197.06 


73.09 


398.63 


0.69 


2092 


29348 


75 


239.19 


28.65 


228.83 


119.40 


418.12 


0.69 


— 


29350 


55 


238.13 


27.69 


242.59 


92.61 


334.90 


0.64 


— 


29587 


207 


239.52 


27.32 


264.53 


365.47 


740.31 


0.87 


2142 


29744 


53 


219.29 


24.72 


261.87 


96.53 


396.59 


0.76 


1939 


30391 


68 


248.60 


26.75 


209.59 


93.25 


271.28 


0.73 


— 


32006 


71 


174.33 


55.04 


170.64 


72.95 


420.65 


0.53 


1396, 1400 


32663 


51 


225.59 


21.36 


186.20 


61.37 


435.82 


0.73 




32909 


79 


226.20 


28.44 


173.18 


69.08 


560.03 


0.56 


2022 


32976 


64 


207.90 


9.50 


195.28 


72.04 


510.21 


0.68 


1808 


33082 


77 


177.27 


54.72 


178.28 


79.96 


382.43 


0.52 


1383, 1396, 1400 


33739 


79 


214.39 


8.24 


172.84 


65.40 


517.22 


0.49 


1890 


33851 


138 


230.43 


20.75 


121.03 


74.65 


354.61 


0.77 




34513 


53 


225.86 


7.88 


261.94 


94.18 


426.82 


0.75 


2020 



Notes. Columns are as follows: 1: ID of the cluster; 2: the number of galaxies in the cluster, A'gai; 3^: cluster right ascension and declination; 5: 
cluster comoving distance; 6: cluster total luminosity; 7: rms velocity of galaxies in the cluster; 8: cluster virial radius; 9: Abell ID of the cluster. 
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Table 1. ... continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


ID 




R.A. 


Dec. 


Dist. 




O" 




Abell ID 






[deg] 


[deg] 


[/!-' Mpc] 


[lo'^r^Lo] 




[r' Mpc] 




34726 


145 


228.09 


7.44 


135.44 


121.20 


506.15 


0.74 


2028, 2033, 2040 


34727 


256 


227.74 


5.78 


234.02 


351.58 


825.94 


1.25 


2028, 2029, 2033, 2040 


35037 


79 


236.16 


36.16 


195.74 


76.68 


691.85 


0.58 


2122, 2124 


36861 


66 


192.32 


55.05 


248.22 


110.73 


493.86 


0.84 


1616 


38087 


169 


168.12 


40.60 


224.46 


209.25 


541.59 


0.84 


1173, 1187, 1190, 1203 


39489 


166 


230.70 


27.80 


212.87 


188.19 


1061.85 


0.72 


2056, 2065 


39752 


108 


231.90 


28.72 


195.69 


116.96 


514.78 


0.75 


2073, 2079 


39914 


63 


233.12 


28.02 


218.38 


69.77 


446.50 


0.65 


2089 


40520 


52 


188.89 


16.52 


209.16 


56.00 


486.67 


0.67 


1569 


40870 


118 


190.34 


18.56 


215.43 


153.64 


717.36 


0.69 


— 


42481 


57 


246.77 


14.17 


151.19 


55.26 


354.04 


0.53 


— 


43336 


68 


207.98 


46.33 


186.81 


73.08 


471.40 


0.60 


— 


43545 


51 


118.25 


29.41 


179.41 


42.75 


577.85 


0.51 


602 


43966 


74 


196.06 


19.31 


190.82 


76.03 


613.87 


0.52 


1668, 1669 


44471 


113 


169.12 


29.26 


141.53 


90.86 


464.08 


0.55 


1185, 1213 


47492 


74 


146.77 


54.50 


138.69 


57.06 


458.67 


0.56 


— 


48448 


55 


140.25 


54.85 


136.34 


41.93 


340.67 


0.54 


784 


50129 


52 


205.43 


26.39 


225.02 


61.83 


449.92 


0.51 


1775 


50631 


101 


116.97 


18.60 


142.30 


86.34 


636.46 


0.55 


— 


50647 


52 


232.32 


52.88 


217.84 


62.30 


525.70 


0.52 


— 


50657 


55 


116.50 


18.24 


153.75 


49.96 


555.58 


0.42 


— 


52913 


67 


129.94 


28.93 


237.98 


116.60 


368.57 


0.77 


690 


56571 


55 


223.57 


54.27 


294.12 


139.89 


457.04 


0.60 


1999, 2000 


57317 


118 


216.74 


16.76 


158.70 


105.37 


516.37 


0.59 


1913 


58101 


105 


168.64 


2.55 


228.30 


122.05 


614.45 


0.92 


1205 


58305 


167 


223.28 


16.76 


135.12 


120.04 


401.84 


0.61 


1983 


58323 


64 


150.66 


32.72 


153.13 


37.74 


393.28 


0.40 


— 


58604 


58 


170.39 


2.86 


150.50 


39.84 


528.86 


0.42 


— 


59794 


90 


187.45 


11.71 


255.93 


141.55 


650.79 


0.75 


1552 


60539 


107 


184.40 


3.65 


230.10 


136.74 


830.47 


0.55 


1516 


61613 


77 


175.58 


7.85 


257.83 


138.13 


517.01 


1.06 


1358 


62138 


124 


223.64 


18.67 


175.90 


127.48 


456.42 


0.78 


1991 


63361 


72 


205.59 


2.26 


230.57 


103.06 


646.10 


0.64 


1773 


63757 


87 


186.99 


8.79 


267.64 


155.17 


653.32 


0.65 


1541 


63949 


80 


179.25 


5.06 


226.17 


112.03 


661.43 


0.55 


1424 


64635 


109 


214.18 


1.95 


162.93 


96.38 


489.43 


0.73 


— 


64702 


64 


184.65 


5.20 


227.92 


84.17 


539.62 


0.74 


1516 


67116 


80 


208.29 


5.17 


236.99 


114.14 


651.61 


0.44 


1809 


67297 


95 


127.24 


30.48 


150.20 


86.98 


770.39 


0.45 


671 


68376 


106 


230.35 


30.63 


230.18 


157.52 


671.54 


0.51 


2061, 2067 


68625 


92 


231.04 


29.90 


336.50 


301.28 


874.84 


0.79 


2069 


73088 


141 


215.50 


48.40 


213.15 


184.99 


631.93 


0.71 


1904 


73420 


68 


122.46 


35.12 


247.01 


105.76 


555.48 


0.73 


628 


74783 


65 


155.30 


23.96 


120.16 


46.02 


401.86 


0.48 





Notes. Columns are as follows: 1: ID of the cluster; 2: the number of galaxies in the cluster, A'gai; 3^: cluster right ascension and decUnation; 5: 
cluster comoving distance; 6: cluster total luminosity; 7: rms velocity of galaxies in the cluster; 8: cluster virial radius; 9: Abell ID of the cluster. 
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Table 2. Results of the tests 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


ID 


Vpec 


Vpecj- 


N30 


un3D 


Pa 


Pa 


N2D 


unlD 


Pe 


Pad 


Psw 


Pkurt 


Pskew 



18 


387.17 


0.753 


2 


1 ■ 10"^ 


4 ■ 10"5 


323 


237.15 


0.857 


3 


0.001 


0.180 


608 


-228.43 


-0.429 


1 


1 ■ IQ-^ 


0.884 


748 


613.95 


0.821 


1 


1 ■ 10"' 


0.007 


793 


-384.68 


-0.747 


3 


0.019 


4 ■ 10"5 


880 


400.67 


0.974 


3 


1 • 10"^ 


0.002 


914 


-704.39 


-1.071 


5 


0.010 


4 • 10"5 


1469 


-342.33 


-0.817 


2 


0.005 


0.198 


1573 


1228.28 


1.651 


3 


0.004 


0.334 


1944 


-122.77 


-0.279 


3 


0.003 


0.020 


2067 


-788.08 


-1.373 


2 


1 ■ 10"' 


0.007 


3714 


-6.37 


-0.018 


2 


0.013 


4 ■ 10~5 


4122 


-1091.85 


-1.133 


3 


0.002 


4 • 10"5 


4713 


-78.38 


-0.123 


3 


0.016 


0.003 


4744 


147.05 


0.298 


4 


1 ■ 10"' 


0.001 


4992 


303.46 


0.489 


2 


1 ■ 10"' 


0.048 


5217 


-19.11 


-0.033 


1 


1 ■ 10"' 


0.216 


7102 


326.16 


0.700 


2 


1 ■ 10"' 


0.022 


7932 


-111.73 


-0.270 


2 


1 ■ 10"' 


0.089 


9029 


-341.66 


-1.095 


2 


0.005 


0.094 


9350 


-129.46 


-0.269 


4 


1 • 10"' 
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0.900 
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0.252 
0.907 
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0.271 
0.850 
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0.851 
0.118 
0.700 
0.430 
0.787 
0.451 
0.221 
0.561 



Notes. Columns are as follows: 1: ID of a cluster; 2: peculiar velocity of the main galaxy (kms^^), 3: normalised peculiar velocity of the main 
galaxy, Vpec/<Tv; 4: the number of components in 3D in a cluster, N^omp', 5: the median uncertainty of classification; 6: /?- value of DS test; 7: p-value 
for a test; 8: the number of components in 2D in a cluster, Ncomp', 9: the median uncertainty of classification in 2D; 10: p-value of /3 test; 11: 
p-value for AD test; 12: p-value of SW test; 13: p-value for kurtosis test; 14: p-value for skewness test. 
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Notes. Columns are as follows: 1: ID of a cluster; 2: peculiar velocity of the main galaxy (fcm.v"'), 3: normalised peculiar velocity of the main 
galaxy, Vpec/cTy; 4: the number of components in 3D in a cluster, A'comp; 5: the median uncertainty of classification; 6: value of DS test; 7: p- value 
for a test; 8: the number of components in 2D in a cluster, Ncomp', 9: the median uncertainty of classification in 2D; 10: p-value of /3 test; 11: 
p-value for AD test; 12: p-value of SW test; 13: p-value for kurtosis test; 14: p-value for skewness test. 



